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Structural transitions in vertically and horizontally coupled parabolic channels of 

Wigner crystals 
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Department of Physics, University of Antwerp, Groenenborgerlaan 171, B-2020, Antwerpen, Belgium 

Structural phase transitions in two vertically or horizontally coupled channels of strongly interact- 
ing particles are investigated. The particles are free to move in the i-direction but are confined by a 
parabolic potential in the y-direction. They interact with each other through a screened power-law 



potential (r n e 



-r/X\ 



In vertically coupled systems the channels are stacked above each other in the 



direction perpendicular to the (x, y)-plane, while in horizontally coupled systems both channels are 
aligned in the confinement direction. Using Monte Carlo (MC) simulations we obtain the ground 
state configurations and the structural transitions as a function of the linear particle density and 
the separation between the channels. At zero temperature the vertically coupled system exhibits a 
rich phase diagram with continuous and discontinuous transitions. On the other hand the vertically 
coupled system exhibits only a very limited number of phase transitions due to its symmetry. Fur- 
ther we calculated the normal modes for the Wigner crystals in both cases. From MC simulations 
we found that in the case of vertically coupled systems the zigzag transition is only possible for low 
densities. A Ginzburg-Landau theory for the zigzag transition is presented, which predicts correctly 
the behavior of this transition from which we interpret the structural phase transition of the Wigner 
crystal through the reduction of the Brillouin zone. 



PACS numbers: 05.20.-y, 64.60.F-, 81.30.-t 



I. INTRODUCTION 

Self-organized systems are of fundamental importance 
in different areas of physics, can be seen in the acceler- 
ated progress in that area. It all started in 1934 when 
Wigner— surmised that, if considering an electron gas, 
where the electrons have no kinetic energy, as occurs in 
low density systems, then these electrons "would settle 
in configurations which correspond to the absolute min- 
ima of the potential energy. These are closed-packed lat- 
tice configurations, with energies very near to that of the 
body-centered lattice. " This was the first prediction about 
self-organized systems of electrons and charged particles, 
which today is known as Wigner crystals. 

At the present it is known that Wigner crystals 
form a body-centered cubic (BCC) lattice in three- 
dimensional^ (3D) space, a triangular lattice in two- 
dimensional^ (2D) systems, while in one-dimensional 
(ID) systems, the energetically more favorable array is 
given by an evenly spaced lattice. Wigner crystals have 
been the study object of several experiments^— and the- 
oretical worksii~— in recent years. 

Recently, the transition between ID and 2D sys- 
tems has been analyzed, known as quasi-one-dimcnsional 
(Q1D) systems. Several studie o 11 ' 17 have investigated 
different properties of the Wigner crystal in that regime. 
When increasing the density, the linear chain structure 
undergoes a zigzag transition that occurs always through 
a continuous transition. This zigzag transition has been 
observed in different cxperiment a 18 ' 19 and analyzed in 
different systems: colloidal systems^—, vortices in su- 
perconductor/ferromagnetic hybrids^. They have also 
been proposed as a possible implementation of a quan- 
tum information processor—. 

The zigzag transition for a system of Coulomb parti- 
cles subjected to a parabolic confinement potential^ was 



analyzed through a Ginzburg-Landau theory. Later this 
theory was applied to more complex particle interactions 
and generalized for a power-law confining potential^. In 
the latter work it was shown that this theory leads to 
a better representation of the phase transition behavior 
close to the critical point. 

In the present work we analyze the phase transitions 
between two vertically and horizontally coupled Q1D 
channels as function of the inter-channel distance and lin- 
ear particle density in each channel. Further the zigzag 
transition of the vertically coupled channels is compared 
with the analytical results obtained from the Ginzburg- 
Landau approach. 

The present paper is organized as follows. In Sect. |TT] 
the model system is formulated. In Sect. IHII the numer- 
ical results are presented with an analysis about the be- 
havior of the system in the different structural phases. In 
Sect. HVI the Ginzburg-Landau Lagrangian is derived and 
the theoretical results are discussed and compared with 
our numerical simulations. Our conclusions are given in 
Sect. H 



II. MODEL SYSTEM 

In this work we consider a system consisting of two 
quasi-one-dimensional channels each with N identical 
particles with mass m and charge q, which are restricted 
to move in the (x, y)-plane. The system is subjected to 
an external magnetic field B, which is applied in the 
direction perpendicular to the channels. The particles 
are confined by a one-dimensional potential limiting their 
motion in the y-direction in each channel. The channels 
are separated by a distance d y from each other in the 
y-direction and a distance d z in the z-direction as shown 
in Fig. HJ 
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FIG. 1. (Color online) Schematic representation of the model 
system, where the filled circles represent the particles in each 
channel, d y and d z indicate the distance between the channels 
in the y- and z- axis respectively. The darker colors in each 
plane represent the increases of the strength of the confine- 
ment potential in the channel. 



The charged particles interact through a repulsive in- 
teraction potential. The kinetic energy of the system is 
given by: 
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where pki = mv^ + qAki represents the canonical mo- 
mentum of particle i in channel k, which is moving with 
velocity Vki- Aki is the vector potential related to the 
magnetic field through the relation B = V x Aki- The 
total energy of the system is given by: 
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where the confinement and the interaction potential are, 
respectively, given by: 
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Here, VJ, a4r (r) represents the inter-particle potential 
which is taken as a screened power-law potential, which 
will allow for the simulations of both long- and short- 
range interactions, as follows: 
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In general rAi.Bj = \*Ai — ?Bj \ represents the relative 
position between the i-th particle in channel A and the 
j-th particle in channel B. The exponent n is an inte- 
ger and e is the dielectric constant of the medium where 
the particles are moving in. In the above, R is an arbi- 
trary length parameter which we introduced to guarantee 
the correct units. The magnetic field is taken to be con- 
stant and is applied in the direction perpendicular to the 
plane formed by the channels, through the vector poten- 
tial defined by = —By^e^. The total energy can be 
written in dimensionless form as follows: 
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where the dimensionless frequency is given by 
v = Vt/uJo, while measures the strength of 
the confinement potential and to = 1/wo i s the 
unit of time. The unit of the energy is given by: 

E = (m Wo 2 /2) n/( " +2) ( g 2 /e) 2/(n+2) i? 2 ("- 1 )/("+ 2 ), 
the distances arc scaled with tq = 
(2g 2 /mw 2 e) 1/(n+2) i?(™- 1 )/(™+ 2 ), and the strengths 
of the magnetic field and the vector potential are 
measured in units of Bq = 2mujo/q and Aq = BqTq 1 
respectively Additionally, the dimensionless parameter 
K = ro/X represents the screening parameter of the 
potential. In order to describe better the behavior of 
the transitions in the system, we define a dimensionless 
linear density 77 as the number of particles per unit of 
length along the unconfined direction in each channel. 
This density is chosen to be equal in both channels. 



III. NUMERICAL MODEL 

In order to understand the behavior of coupled sys- 
tems, we consider two different cases. Case A corre- 
sponds to two vertically coupled channels where the dis- 
tance between the channels in the y-dircction is taken 
zero. In this case the channels are arranged above each 
other. Case B consists of two horizontally coupled chan- 
nels with d z = 0, which are aligned parallel to each 
other separated by a distance d y . In order to find the 
ground state configuration for each case we have per- 
formed Monte Carlo simulations, complemented with the 
Newton optimization technique's. Using those methods 
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FIG. 2. Phase diagram of vertically coupled channels as a 
function of the linear particle density and the distance be- 
tween the channels in the z-direction. Solid and dashed lines 
represent first and second order transitions, respectively. The 
different phases are discussed in the text and shown in Fig. [3j 



we construct the phase diagram of the ground state con- 
figuration of the system at zero temperature as a function 
of the separation between the two channels and the linear 
density. 



A. Case A: Vertical Coupling 

The phase diagram of vertically coupled channels is 
shown in Fig. [2] for an interaction potential with param- 
eters n = 1 and k = 1 and fixed value of the frequency 
v = 1. The phases are represented by the number in 
each delimited region, the solid lines represent first order 
transitions and the dashed lines second order transitions. 

From this phase diagram we see that the zigzag transi- 
tion occurs only between the one-chain configuration and 
phase 1 where particles of each channel are organized in 
a single chain, this transition is continuous and of second 
order. The different phases are shown in Fig. [3] One can 
see that the activation point , defined as the point where 
the one-chain arrangement is no longer the only possible 
ground state configuration for any value of d z , is given 
by r\act = 0.447302 for d z = 0. This value of the linear 
density represents exactly half the value of the critical 
density (j] c ) for the zigzag transition in the case of a sin- 
gle Q1D channel as was found in Refs. |Tl| and[28|. Notice 
that for increasing d z , the critical density increases as 
well, which means a change in the total dimensionality 
of the system. Additionally, one can observe that for 
d z > 0.5, r\ c increases more quickly than for d z < 0.5. 
This implies that the system tends to behave as two de- 
coupled channels, until the point where the zigzag tran- 
sition disappears. 

The different ground state configurations identified in 
Fig. [3] are shown in Fig. [31 In this figure, the red 
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FIG. 3. (Color online) Configurations of the different crys- 
talline phases of the particles as enumerated in the phase dia- 
gram presented in Fig. [2] The red solid (blue open) symbols 
are the particles in the top (bottom) channel. 



filled circles represent the particles in channel A and the 
blue open circles the particles in channel B. The gray ar- 
rows between two configurations indicate that the tran- 
sition between them occurs continuously, following the 
process described through the small black arrows inside 
each phase configuration. The small black arrows above 
each particle represent the displacement direction when 
the separation between the channels is increased. It is 
interesting to observe that the phases 7, 8 and 13 are 
composed of several configurations whose transition oc- 
curs through a continuous phase transition indicated by 
the small black arrows plotted inside the pictures. 

The eigenvectors of the first non-zero eigenfrequency 
for each ground state configuration is shown in Fig. 3) 
The direction of the eigenvector is of particular interest 
because it contains information about the onset of melt- 
ing in that crystal, when temperature is increased. 

In order to know the evolution of the system it is in- 
teresting to show the behavior of the total energy as a 
function of the control parameter. This can be seen in 
Fig. [SJ where the contour plots of the energy per particle 
are plotted as a function of the linear density and the 
vertical separation between channels. The color scale in- 
dicates the magnitude of the energy and the lines show 
the iso-energy lines on the plot. Notice that by increas- 
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FIG. 4. (Color online) Eigenvectors of the first non-zero fre- 
quency of each phase for vertically coupled systems. 
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FIG. 5. (Color online) Contour plots of the energy per particle 
as a function of the linear density r\ in each channel and the 
vertical separation d z between the channels. The red lines 
represent the iso-energy lines of the system. 



ing the density, the iso-energy lines become increasingly 
curved showing the decoupling of the channels, due to 
the increase of the interaction between particles in the 
same channel, while the energy contribution from the in- 



teraction between channels becomes significantly smaller 
with increasing d z . 

With the aim to explain the order of the transitions, 
we plot in Fig.|6]the first non-zero eigcnfrequcncy for dif- 
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FIG. 6. (Color online) (a) First non-zero eigenfrequency; (b) 
Energy per particle; (c) First derivative of the energy; (d) 
Second derivative of the energy. In all figures we have plotted 
as a function of the vertical separation between the channels. 
The colors showed at the top represent the different densities. 
The inset in (d) shows the distance of the chains from the 
a;-axis of the channels as a function of the vertical separation 
for T) = 1. 21. 
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FIG. 7. (Color online) Phase diagram of the system as a 
function of the linear density (rj) and the distance between 
the channels in y-direction (d y ). All phase transitions are 
first order transitions. 



ferent densities. We also plot the energy per particle and 
its two first derivatives with respect to the vertical sepa- 
ration. In all cases the figures show the behavior of these 
quantities as a function of the separation between chan- 
nels, and the value of the density taken for each curve 
is printed at the top of the figure. From the smooth and 
monotonic behavior of the energy in Fig. |6|b) and the 
jumps in its derivatives (see Figs. |6][c,d)), wc can de- 
duce the order of each phase transition, which following 
the Ehrenfest classification^ are given by the order of 
the lowest derivative of the energy which exhibits a dis- 
continuity. From Fig. [H^a) we note that the transitions 
between phase 1 and phase 2 (shown for rj = 0.72 around 
d z = 0.798) and between phase 2 and phase 3 (shown for 
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FIG. 8. (Color online) Contour plots of the energy per parti- 
cle as a function of the linear density rj in each channel and 
the horizontal separation d y between channels. The red lines 
represent the iso-energy lines of the system. 



rj = 0.72 around d z = 0.866) are of first order which can 
be seen by the jumps in the first derivative of the energy 
at the transition points. The transitions of second order 
are characterized by the softening of an eigenmode, i.e. 
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FIG. 9. (Color online) (a) First non-zero eigenfrequency; (b) 
Energy per particle; (c) First derivative of the energy; (d) 
Second derivative of the energy. All figures are plotted as a 
function of the horizontal distance between the channels. The 
colors showed at the top represent the different densities. The 
inset in (d) shows the distance of the chains from the center 
of each channel as a function of the horizontal separation for 
rj = 1.19. 
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u> = 0. The inset in Fig. EJd) shows the y-position of 
each chain for r\ = 1.21, where one can see the continu- 
ous transition from 3 to 6 chains, the first order transition 
between 6 and 4 chains and the second and continuous 
transitions between 4 and 2 chains. The red and blue 
lines represent the positions of the particles in channel A 
and channel B. 



small black arrows have the same meaning as previously. 
Notice that phases 14 and 15 consist of several phases. 
The transition between those phases occurs through con- 
tinuous transitions. In Fig.fTUTb') the normal mode of the 
first non-zero frequency is plotted for each configuration. 
IV. GINZBURG-LANDAU MODEL FOR 
LINEAR-ZIGZAG TRANSITION 



B. Case B: Horizontal Coupling 

In this part we analyze horizontally coupled channels 
(d z = 0) for the same parameters as in case A. The phase 
diagram in this case is shown in Fig. m function of 
the linear density r\ and the vertical separation between 
the channels d y . 

From this phase diagram we notice that the zigzag 
transition docs not occur in horizontally coupled chan- 
nels. The one-chain configuration only exists when d y — 
and for r/ < r) act . Additionally the small number of 
phases is a result of the weaker interaction between chan- 
nels. 

In Fig.[5]we show a contour plot of the energy per parti- 
cle and iso-energy lines as a function of the linear density 
and the horizontal separation between the channels. No- 
tice the smooth and monotonic behavior of the iso-energy 
lines, which align parallel to each other with increasing 
linear density. When increasing d y , the interaction be- 
tween the channels decreases quickly due to the separa- 
tion between channels which is directly related with the 
confinement potential as can be seen in Eq. ([5]) , where d y 
contributes to the total energy of the system introducing 
a linear term of the relative position between particles 
in the y-direction. This contribution leads to the disap- 
pearance of the harmonicity of the system. 

In Fig. IHl^a) we show the first non-zero frequency as 
a function of the horizontal separation for different den- 
sities, which shows clearly the phase transition points. 
The monotonic and continuous behavior of the energy 
(Fig. Elb)) and the jumps in the first derivative of the 
energy, as shown in Fig. |U(c), indicates, that all phase 
transitions are of first order. In those transition points 
we also can see that the first non-zero eigcnfrcqucncy has 
a jump, which shows the abrupt change in the vibration 
mode of the system in the phase transitions. As an illus- 
tration about the first order transitions in the system, we 
show in the inset of Fig. Od) the y-distance of the chains 
from the center of each channel, where one can see the 
transitions between 3, 8 and 4 chains for r\ = 1.19. In 
this figure the red and blue lines represent the particles 
in channel A and channel B, respectively. 

The different configurations of the phases shown in 
the phase diagram of the horizontally coupled system, 
are plotted in Fig. [TUTa - ). where the gray arrows and the 



The results of previous section indicates that the tran- 
sition between the one-chain arrangement and the zigzag 
configuration occurs in the low density region and only in 
the case of a vertically coupled system. To analyze this 
transition (i.e. d y = 0), we start by considering the sys- 
tem in the situation that the one-chain configuration is 
stable in each channel and close to the zigzag transition 
point. The equilibrium positions of all particles are given 
by rf™ = (x«» 0, (k - 1)4) with a;«» = [2i + (2 - k)]/2 V , 
where k represents the channel number (i.e. 1 or 2) and 
i the particle number in that channel. Further we con- 
sider small oscillations around the equilibrium positions 
in the y-direction. The position of each particle becomes 
thus Tki = rj™ + Ar^i , where the small displacements are 
given by Ar^ = (0,yfci,0). Now, ensuring that the par- 
ticles are oscillating around these ground state positions, 
and assuming that these oscillations are much smaller 
than the distance between the two nearest particles, we 
expand r\i,2j into a power series. 

As found numerically, this transition is continuous and 
the configuration after the transition point consists of a 
zigzag configuration (phase 1 in Fig. [2]), but with all 
particles forming a single off-center chain in each chan- 
nel (i.e. yii = c and y2i = — c, with c as a real number). 
From that configuration we can observe that there is a 
competition between the interaction potentials Vii, V22 
and V12 which are defined as the last three terms in Eq. 
(|5j) respectively. The subindexes represent the interact- 
ing channels. The first two interaction potentials (Vn, 
V22) do not contribute to the solution in this case, since 
all particles in the same channel are aligned parallel to 
the x-axis. Thus, the interaction between the particles 
in the different channels, as well as the confinement po- 
tential, are responsible for this transition. 

Following Ref. H^, we obtain the Lagrangian density 
of the system in the continuum limit by considering only 
the interaction potential V12: 
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h 2 (n, k, c z ){d x ip(x)) 2 + 5 v (n, k, c z )ip 2 (x) 
+A(n,k,c z )ip 4 (x) , (6) 



where ft = k/ij, c z = r]d z , S v (n,k 7 c z ) = 2{v 2 + B 2 /4) — 
w 2 {n, h, c z ) and the order parameter ip{x) represents the 
separation between channels (i.e. ip — 2c). The coeffi- 
cients arc given by: 
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FIG. 10. (Color online) (a) Configurations of the crystalline organization of particles enumerated in the phase diagram presented 
in Fig. [7] (b) Eigenvector of the first non-zero frequency of each phase configuration. 
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where i?j = ^/ (j + 1/2) 2 + c 2 and fco = fco/f7 represent 
the value of the wavevector in the stability point. 



A. Stability Point 

The stability point is calculated numerically, as the lo- 
cation of the minimum of the second order term of the 
total energy, and the value of the wavevector is called 
fco- Then the stability point is found from the condition 
min<J w (n, re, c z ). Notice that the value of fco depends on 
the distance between channels. For c z = we can find 
analytically that fco = 27r, which implies that the particle 
density (77) of the system in each channel will be reduced 
to half of the density required to find the stability point 
in one channel as found previousl y 11 ' 28 . Next, from the 
stability condition we find for v = 1 that the activation 
point (minimum value of density to find the zigzag tran- 
sition) is 7? c ( mi „) = 0.447302. 

We show in Fig. QT] the behavior of fco as a function of 
the particle density in each channel. Notice that for den- 
sities lower than r] c r m in)i which represents the region of 
the one-chain configuration, the value is fco = 2-7T. This 
value indicates that the first Brillouin zone is twice as 
large as in the one channel case. This value is related 
to 77, because the total linear density of the system is 2rj 
when c z = and it decreases by increasing rj. For lower 
densities {rj < r) c (min)) the linear arrangement is pre- 
served due to the weaker interaction potential between 
the particles as compared to the confinement potential 



and this for any separation between the channels. How- 
ever, when rj > ?7 e ( TO j ra ), the zigzag configuration is found 
to result in a lower interaction energy between the par- 
ticles in the different channels. This behavior can be 
reduced by increasing c z and therefore reducing by the 
total density of the system. In that situation fco is of ma- 
jor importance, because while the separation increases 
and the total density decreases the first Brillouin zone 
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FIG. 11. Wavevector corresponding to the second lowest 
eigenvalue as a function of the particle density in each chan- 
nel. System parameters: v = 1, n = 1, and k = 1 
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FIG. 12. (Color online) Critical value of distance between 
channels as a function of the particle density in each channel. 
The solid line is the theoretical result and the open circles are 
the results from our MC simulations. We took the parameters 
v = 1, n — 1, and k = 1 



FIG. 13. (Color online) Distance between channels as a func- 
tion of the vertical separation between the channels, where 
the linear density in each channel is r\ = 0.52. We took the 
parameters v = 1, n = 1, K = 1. 



becomes smaller reducing the dimensionality of the crys- 
tal in each channel. This behavior is shown in Fig. QT] 
where one can see the monotonic decrease of ko with 77, 
when the linear density exceeds the critical density. 

In order to find the ground state configuration, we 
minimize Eq. and we obtain the equation of mo- 
tion of the system. Then, after minimization, the so- 
lution is reduced to the well-known analytical expression 
tp = \J—8 v (n, k, c z )/2A(n, k, c z ) which is valid when the 
density is lower than the critical density, as was found in 
Ref. [2S| for parabolic confinement. 

Fig.[l2]rcprcsents the critical separation between chan- 
nels as a function of the particle density. This critical 
value corresponds to the separation c z where the zigzag 
transition takes place. This critical separation is calcu- 
lated as the lower value of c z at which -0 = for a given 
rj. For points below the curve in Fig. [12] the zigzag ar- 
rangement is the ground state configuration. The open 
red circles are the results for the zigzag transition points 
calculated from our MC simulations. The plotted region 
corresponds to the region where the zigzag transition is 
allowed, as was analyzed in previous section. One can see 
that there is perfect agreement between our theoretical 
results and the results from the MC simulations. 

In Fig. [13] we show the distance between chains in the 
confinement direction as a function of the separation be- 
tween the channels c z for a linear density 7/ = 0.52. The 
open circles show the results from our MC simulations 
and the solid line represents the value of the order pa- 
rameter tp calculated from the present theory. The good 
agreement indicates that the Ginzburg-Landau theory is 
able to describe correctly the behavior of the system close 
to the zigzag transition point. 



V. CONCLUSIONS 

We studied two coupled Q1D channels, consisting of 
a system of interacting charged particles confined by a 
parabolic trap in each channel, where the linear particle 
density in both channels is given by 77. The ground state 
configuration at zero temperature was analyzed. 

The structural transitions between phases were studied 
as a function of the linear density and the separation be- 
tween the channels. We found a very rich phase diagram 
in case of vertically coupled channels, with first and sec- 
ond order transitions. The horizontally coupled system, 
on the other hand, exhibits a very restricted number of 
phases and all transitions are of first order. The latter 
can be traced back to the linear term in the relative dis- 
tance between particles in the different channels, which 
appears in the energy expression. This linear term results 
in a rapid decrease of the interaction between channels. 

Our simulations show that the zigzag transition occurs 
only in case of vertically coupled channels, which repre- 
sent another effect of the linear term in the energy. In 
order to understand the behavior of this zigzag transition, 
we derived a Ginzburg-Landau equation and determined 
the behavior of the system close to the zigzag transition 
point. We analyzed the order parameter and its depen- 
dence on the linear density and the separation between 
channels. From theory we found that, for vertically cou- 
pled channels, an increase of the density produces a re- 
duction of the first Brillouin zone of the Wigner crystal 
after the activation point rj act . Within the density range 
where the zigzag transition is allowed, the vertical sepa- 
ration could be used as a tunable parameter to modulate 
the transition point. This result was found from theory 
and there is a perfect agreement with our numerical re- 
sults, which shows that the presented theory is sufficient 
to understand the behavior and the nature of the zigzag 
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transition in Q1D coupled channels. 
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